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SECTION I 


NEAR FIELD TO FAR FIELD TRANSFORMATIONS 


1.0 Introduction 


1.1 Background 

Near field measurements to determine far field patterns of 

antennas are felt to be an effective alternative to far 

12 3 

field range measurements. ' ' 





FIGURE 1 

FAR FIELD ANTENNA PATTERNS FROM 
NEAR FIELD MEASUREMENTS 
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This is especially true for microwave antennas and most 
especially for phased array antennas. Near field measure- 
ments allow measurements to be made indoors. An indoor 
range provides the all-weather capability to enable tests 
to be performed in a controlled environment. In addition, 
large-antenna logistics and mounting problems are simpli- 
fied by near field measurement systems. These systems are 
felt to be time and cost effective, with accuracy compar- 
able to that of far field ranges. The requirements for a 
more complicated and expensive measurement system with 
more expensive probe calibration are not felt by pro- 
ponents to be serious disadvantages. 

1.2 Project Objectives 

As a part of this contract,! the feasibility , of acoustic 
verification of microwave near field to far field trans- 
formation algorithms using the Phased Array. Sonic Simulation 
System was studied. Existing electromagnetic near field 
measurement techniques and transformation algorithms 
(equations) were investigated .''^"'"^ ' It was analytically 

determined that acoustic verification is valid^ We—weuid 
— Ti'ke—to— exp'1'd-e-i-tly point out that .although— farr~~f±eid 

pa-t-te-r-ns— f-rom— near_f:ie-l-d— patterns— for— acoustic— arrays 
6 7 

are— vsrlxd — (arrd"~fn~"f arct~ma-y— ha-ve-be. en the e a r 1 ije,s,t-.us e 

o f— n e a-r— f-i e 1 d— to —far ~ f ±e±d— t r ari s’f ormatf i oiTsTr~ what— we— are 

i'-- 

-t-alkiiig— ab<aut--here— i-s j acoustic simulation of electromagnetic 
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near field to far field transformations A - — fOne merely needs 
to separately simulate acoustically the transmission of the 
two orthogonal polarization components of the EM wave and 
then superimpose the results to get the analytical descrip- 
tion of the total field. } — " 


The acoustic simulation of electromagnetic near field to 

far field transformation is verified for the "modal 

- ! ■ - 4 — 9 — 

expansion" method. ' ' In the modal expansion method, 
data from antenna near field measurements are converted 
to a summation'- or spectrum of modes corresponding to wave 
numbers in the measurement coordinate systems ^tth~'a~'~' 
rec t a n g.u lar— coordinate 's vs~t’emT~'the field at all poTnts~Th 
front Q f„a— plane-~aperture""of"any~~ap'e'rture"drs'tri , butuon 
may— b e— r e g a-r de d—as— arising-f rom^~ali“a'ggre'ga'te-“'Of--pi-a n e 
^aves-~trave-ld~ng~iTf'^ari~gu's~"dire‘ctions-r'" > The —amp-l-i-t-ude^ 
'anU _ pKase' 6*f"th^^ave^'“as _ ^'Tunctidh“o'f~th'e'rr“direct-ion^ 


of^trav-e-1— cons-titu-tes— an-^angular-spectrum.- — ^^Vno-thex— way 
o f—iooking~at — th'e'"same'"'th-ing,— and— thi s— iS‘-the--phx-los op i cal 
a pproac h— taken— in— the-|-veri-f ication me tho^\u-sed— in— -th-i-s- 
repox4-p^is that all the far field pattern information is 
contained in the near field and thus is contained in 
measurements made in the near field. Fourier transformation 
of those measurements preserves the far field information 
is a spectral form that is then readily extractableT~] 


'F'C\4£~ 

A Tt o Vi NX ft ^ <>'\ PavNJ ? ^SJi Ki v 
Vl «\v)UR"' t (J\ 
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1.3 Near Field/Far Field Definitions 


Before presenting the details of the mathematical verifi- 
cation, it is important to define what we mean by near 

g 

field and far field. It is commonly accepted that the 
far field is that region which is at a distance from the 
antenna greater than 


where D is the largest dimension of the antenna and A is 
the wavelength. The near field is the near radiating 
field and is the region between the far field and the 
near reactive field. The near field is considered by 
many to begin at a distance from the antenna greater 
than 

X/2tt . 

However 10 this defines the point where the reactive 
near field equals the radiating near field for a small 
dipole. Suffice to say that for most antennas, that 
point is usually less than a few wavelengths from the 
antenna. Thus near field measurements are usually made 
within a few wavelengths from the antenna. Strictly 
speaking, there is nothing in the verfication proof to 
be developed below that restricts the measurement plane 
to any region or distance. 
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2.0 Near Field to Far Field Transformation - Mathematical 
Verification 


2.1 Wave Equation Basis 


Consider the well known wave equation governing electro- 
magnetic radiation in a source-free region and given in 
terms of the time-varying electric field vector, E: 


V 2 E = ye - 


where V is the Laplacian operator. For a monochromatic 
(ie. , single frequency) field. 


E = Real {Ee 3 } 


where the "hat” denotes the phasor. Let the coordinate 
systems of interest be defined by Figure 2. Substitution 



FIGURE 2 


COORDINATES DIAGRAM 
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of (3) into (2) then yields 

- /"V 

V 2 E = -yew E. (4) 

2 

Let us express E and V in rectangular coordinates: 

A A 

E = E (x,y , z) 



The propagation velocity, c, is given by 


1 



Since 

c = fX, 

where f is the frequency, the wave equation as described 
by (4) can be rewritten as 

v E = -K E (5) 

where k is the propagation constant (wave number) expressed 
by : 

k = 2 tt/A. 

Let us now take a two-dimensional Fourier transform with 
respect to the x and y dimensions of the wave equation. 
Taking such a Fourier transform of both sides of (5) yields 
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( 6 ) 


L 


~ ~ -jk x -jk y 

[vE(x,y,z)Je x dx}e Y dy 


r°° f 00 9 ~ -jk x -jk y 

= . { [-k E (x,y , z) ] e X dx}e Y dy 

J —00 J -*oo 


where k x and k^ can be thought of as spatial frequencies 

Using the rectangular coordinates' expression for the 

2 

Laplacian operator, V , (6) becomes 


C«> ,°° „2 ~ -jk x -jk y 

{ [ — ^ E(x,y,z)]e x dx}e Y dy 

— 00 J —00 J 


-00 J -00 9x 


+ { 


00 r°° 9 2 - _ -jk r y -jk v x 


f 

_ 


E(x,y,z)]e y dy)e x dx 

-» 3y 


f°° r ci 2 


— ooJ 


^ a - j (k x ' + k y) 

[ — j E(x,y,z)]e X Y dxdy 

-oo 3^ 


= -k‘ 


r°° 


E (x,y , z) e 


-j (k x + k y) 


Y dxdy . 


(7) 


Let the double-integral term on the right side of the 
equal sign in (7) be defined as the phasor/vector field 

A 

B: 


~ a r f 

- (k x ,k y ,z) - J_ 


E (x,y , z) e 


-j (k x + k y) 


x 


y dxdy. (8) 


Note therefore that the third term to the left of the 
equal sign in (7) then becomes described by (9): 


- 9 - 


[— j E (x,y , z)] e 

-°° 3 Z 


-j(k x + k v) 


r°° 


3 z‘ 


r a 


E (x,y , z) e 


J dxdy 
~j Ck x + k y) 


dxdy 


— -CQ J — CO 


3 Z 


2 §(k x ,k y , z ) . 


( 9 ) 


Before examining the first and second terms to the left 
of the equal sign in (7), let us look at the Fourier 
transform of differentials. . If the Fourier transform 
of a function, f(t), is given by 


S*{f (t) } 


,°0 

f(t)e ^ wt dt =F(jw), 

— oo 


( 10 ) 


then f(t) is given by the inverse Fourier transform 


f (t) = F (joi) } = ■— 


F ( jw) e^^dw . (11) 


a n 


Therefore f(t)} is obtained as a direct consequence 


dt n 

of taking the n ' th derivative of f(t) with respect to t in 
(11) followed by application of the rule given by (10) : 

,n 

— — f(t)} = ( jw) n F ( ju) = (j<jj) n ^{f (t) }. 
dt 


Thus in the first term of (7) we get 

r 
k. 

dx 


r 00 „2 ~ -jk x 

I r 3 .. J x 


[ — 2 - (X/y» z)] 1 

-°°3x 


7 r 00 ~ -jk X 

= (j k x ) E (x,y,z)e x dx. 

J —OO 
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After taking the Fourier transform with respect to y 
of this result and using (8), the first term of (7) 
becomes 


r~ c 00 g2 


{ [i-y i( X ,y, 2 )]e' :|kxX dx}e" : ’ k y y dy 
_coJ -oo g x 


- 


a 


E (x,y , z) e 


-j (k x x + k y y) 


dxdy 


= -k B (k ,k , z) . 
x — x' y' 


( 12 ) 


In a like manner, for the second term in (7) , we have: 


r°° r°° 2 

j { [— y E(x,y ,z)] e 

j — OO — OO Ti\7 


-jk y 'D k v x 

Y dy } e dx 


-°° 3y 


= - k B (k ,k ,z). 
y — x y 


(13) 


Substituting the results given by (9) , (12) , and (13) 

into (7) gives us the Fourier transformed wave equation 

✓\ 

in terms of B (k ,k , z) : 

— x y 

n A « A -2 A A A 

-k B -k b + B = -k B. (14) 

x y dz Z 

Rearranging terms gives : 

~ B(k,k ,z) = -k 2 B -(k.,k ,z) (15) 

~ x Y Y 

where 
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2 


( 16 ) 


k=k 2 -k-k . 
z z y 


Equation (15) has a solution given by 


B (k x ,k y# z) = A (k x ,k y )e 


jk z 
J z 


(17) 


(Proof is by substitution, which can almost be seen by 
inspection.) (Note that: 


A 


A 


A (k ,k ) = B (k ,k , 0).) 
xy x y 

A 

Since B (k ,k ,z) as given by (8) is seen to be a two- 
— x y 

A 

dimensional Fourier transform of E with respect to x 
and y, then the uniqueness property of Fourier transform 

/v 

pairs means that E(x,y,z) is given by the inverse Fourier 
transform: 


E (x,y , z) 




j(k x + 

(k x ,ky,z )e 


k Y) 

dk x d y 


(18) 


Using the wave equation solution given in terms of B 
expressed by (17), (18) becomes: 


E (x,y , z) = — ~ 

4tt 


r ~ i k 7 z j(k x+k y) 

[a (k ,k ) e z Je x y dk„dk... (19a) 
— x y 


x y 


The result given by (19a) can be expressed in compact 
vector notation as: 


E (r) 




A (k x , 


k )e^ 

y 


£. dk dk 
x y 


(19b) 
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where k=ka +ka + k a , 
— x— x y— y z— z 


r = xa + ya + za 

-x -y -z, - 

and a , a , a are unit direction vectors, 
x “ ’Y z 

2.2 Physical Concepts/Terminology 

Let us now look at the physical significance of the terms 
in the expression for the electric field given by (19b) : 

A 

• E(r) is the electric field, expressed as a vector 
intensity, defined by r, the magnitude and direction 
of that point from origin of the coordinate system. 

A "i k * it 

• A(k x ,ky)e— — is an expression for a plane wave. 

✓\ 

• A(k x ,ky) is a plane wave spectrum . This arises from 
the fact that the form of (17b) shows that E is the 

/\ A 

. inverse Fourier transform of A. Therefore A is the 

A 

Fourier transform (and therefore the spectrum) of E: 


A(k ,k ] 
— x y 


E(x,y,z)e — dxdy, 


00 -00 



(sum) 


A (k ,k ) e-^— — dk dk 

— xv x y 

of plane waves of every 


is a superposition 

"mode" k , k . 

x y 


The field E(r) is therefore expressed as a "modal expansion" 
(a "plane wave expansion") by (19b.) This equation shows 
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that the field may be regarded as being made up of a 
superposition of plane waves. 

2.3 Propagating/Evanescent Plane Waves 

From (19a) and (.19b) we use the relationship 

ik • r j (k x x + K Y \ jk z z 

— ( k x ’ k y ) e ~ “ = [^ (k x ,k y )e Y l e 

We see that the exponential term in z can be expressed as 


jk z 

J z - Yz 

e = e 


where Y is the propagation constant and where from (16) 
we have 



( 20 ) 


In general, Y is composed of real and imaginary parts, 
the attenuation and phase constants, respectively: 


-jk z = Y = ex + jS, 


( 21 ) 


where a and 3 are real numbers. In (16) then, when 


k 2 + k 2 > k 2 , 
x y ' 


k is imaginary. This means that k can be represented as 
z z 


k 2 ° 5 “' 


and therefore 
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e z = e ” az , for k z imaginary. 

Thus we see that those waves for whick k z is imaginary 
will "evanesce" (attenuate rapidly to zero in a direction 
toward the far field). , When 




k z is real and can be represented (choosing the minus value 
for the square root radical) as 



Therefore in this case, 

jk z -j(3z 
e = e . 

n B z 

The term e J is merely a phase term. Thus those waves 
for which k z is real will propagate (will exist at infinite 
values of z, the far field) . 

Thus from (19b) we see that the electric field is made up 
of a superposition of propagating and evanescent plane 
waves . 


2.4 Near Field Measurements 

In a plane parallel to and located at a distance, z q , 

A 

from the aperture plane, the electric field, E as 
represented by (19a) , is given by 
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1 


oo 


oo 


E (x,y , z ) = 


4 TT J — oo J — oo 


~ jk Z Z 0, j (k x X+k V Y) 

[A (k , k ) e °]e y dk__dk_ . 

— x y 


x y 


Therefore the uniqueness of Fourier transform pairs 
dictates that: 


A(k x ,k y ) e 


jk z 
J z o 


■r j 


r - 


E(x,y,z Q )e 


-j (k x + k v) 


y dxdy. (22) 


Thus the plane wave spectrum, A, can be calculated (using 

/\ 

(22)) from measurements of the electric field, E(x,y,z Q ) , 
made in that plane: 


A (k x , k y )= e 


- jk z 
J z o 


- - j (k x + k y) 

E ( x , y , z q ) e y dxdy . (23) 


Equation (23) is valid irrespective of the magnitude of z q . 

A 

Thus A can be determined from measurements made at a distance, 
z q , which puts that plane in the near field of the antenna 
(see Figure 3) . Once we have calculated A from such near 



(can be arbitrarily located 
with respect to antenna) 



FIGURE 3 

ORIENTATION OF MEASUREMENT PLANE 
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field measurements, we can then determine the electric 

field, E, for any r(x,y,z) in Equation (19b). This is, 

✓\ 

we can determine E anywhere by substitution of the result 
of the calculation of A from (23) into (19b): 


E (x,y , z) = 


4 7T ' 


A * V 

A (k , k )e 3 - - dk dk . 

_oo- x y x y 


(19b) 


"Anywhere", of course, includes the far field which is 
what we're especially interested in here. 

A 

The electric field, E(.x,y,z), in the far field (the "far 

A 

field" pattern) can be evaluated by evaluating A(k ,k ) 

x y 

in (19b) at particular values of k^ and k y . The method 

7 

of "stationary phase" can be used to derive the result. 
We will use that method below. But first, let us derive 
some useful relationships. 


2.5 Field Relationships 

A 

The electric field, E, and its corresponding plane wave 

A 

spectrum. A, can be expressed vectorially in rectangular 
coordinates in terms of the unit direction vectors and 
directional phasor components to have (19b) yield: 


E (x,y , z) = a E (x,y , z) + a E (x,y,z) + a E (x,y,z) 
** x y y 





+ 


a 2 A z (x,y,z)]e j ^ 


— dk dk . 
x y 


Therefore we have the triad 


E x (x,y,z) 


4 it 


QQ #00 

~CO J — C 


A„(k ,k )e^— — dk dk (24a) 

oo.-oo x x y X y 


E y (x,y , z) 


4 TT ' 


_A (k *' k y )e3 " " E dk * dk y <24b) 


E z (x,y , z) = — j 
4 IT 


, A z (k x' k y )eJ ' ’ £ dk x dk y (24c) 

— 00 J —00 - 1 - 1 


In a source-free region (which is the case with our near 
field measurements) : 


V • E = 0 

8 8 9 

where V = a„-^— +a- 5 ~ +a . 

—x 9x — ydy — z9z 

Therefore in rectangular coordinates we have 


9E (x,y,z) 9E (x ,y , z) 9E (x,y,z) 
x _ j . Z_ z 


9x 


9y 


9 z 


= 0 (25) 


Substituting into (25) for the components of E in terms 

A 

of those for A as given by (24) yields 


4t t" 


— OO * 


A (k ,k ) „ 
_oo x x y 9x 


r j (k x + k y + k z)-, 

U Y J dk ~ dk 


x y 


+ 


4tt‘ 


_„V k *'Vw [ ' 


j (k x + k r y + k z) 


y z ] dk dk 


X y 


r 00 r* 


4tt‘ 


— OO' — . 00 


WVA [e 


5 r. jl V + V + k * z) iak_.dk 


] "X y 


0, (26 
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which, upon performing the partial dif fenentiations , becomes 


4 7T ‘ 


— CO j —00 


[jk A (k ,k ) + jk A (k ,k ) 
L J x x x y y y x' y 


+ jk A (k ,k )]e 
J z z x y J 


j (k x + k y + k z) 

J v 7 


X 


dk dk = 0 (27) 
x y 


Cancelling out the j’s in (27) and defining the resulting 

A 

expression as f(x,y,z) gives us 


~ i r 03 ^ jk z 

f(x,y,z) = — - [k • A (k x ,k )e z ] 

4 7T J — oo J — oo * 


D (k x + k v y) 

e y dk dk = 0. 

x y 


(28) 


Because of the uniqueness of Fourier transform pairs, 


/v jk z 

k • A (k , k ) e Z = 

x. y 


f (x,y , z) e 


- j (k x + k v) 


Y dxdy . (29) 


But from (28) , 

f (x,y , z) = 0 
Therefore 

^ jk z 

k • A(k ,k )e =0 for all z. (30) 

x y 

jk z z 

But e is in general not equal to zero for all z. 

Therefore in order for (30) to be true for any z, we have 
the result that in a source-free region 

A 

k • A = 0.. (31) 
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This means that 


kA + k A + k A = 0 
xx y y z z 


(32) 


from which we get 


/V -i 

A = - r— [ k A + k A 1 
z k L z x y y J 
z J 


(33) 


where k = /k 2 - k.. 2 - k_ 2 


x 


and where (using (19)) : 


A = 
x 


E^e ^dxdy 


-jk z 


= e 


z o 


E x (x,y,Z Q )e 


- j (k x x + k r y) 


y dxdy (34a) 


A = 

Y ; -co; 


E e — dxdy 

—co y 


= e 


-jk Z r°° coo 
Z O ’ 


E (x ,y , z ) e 

— CO J —CO 


- j (k X + k y) 


y dxdy (34b) 


The key result is that since, from (33) , A is a function 

z 

/V A 

of A x and A^., measurement is only required of the x and y 

/\ /v 

components of E to get A. Stated in more general terms, 

A 

only two measurements of E, each taken at different 
orientations of the measurement antenna, are required to 

A 

obtain all of the components of A. 
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The vector fields, A and E, can also be expressed in 
spherical coordinates as 

A A A A 

A = a r A r + a Q A e + a^ (35a) 

and 

A A A A 

E = a E + a Q E Q + a,E, (35b) 

— — r r 0 6 — tj) 4> 

where a , a„ , a,, are the unit direction vectors. In the 
— r — u — (p 

far field of an antenna the radial component of the electric 

/\ 

field, E r , is negligible compared to the other components. 
Thus in the far field 

E = ^ 0 E 0 + vv (36) 

2.6 Far Field Solution Using Method of Stationary Phase 

The far field solution for the electric field, E, in (19b) 
is obtained below after R. Collin 11 using the method of 
"Stationary Phase". 

The geometry of the situtation is depicted by the coor- 
dinates diagram in Figure 2. In rectangular coordinates, 

r = xa + ya + za . (37) 

— —x — z — z 

It is easy to show from Figure 2 that the x,y, z components 

of r are given in terms of the spherical coordinates, r, 

0 , and 0 as 
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X 


(38a) 


= r sin 0 cos 4> 
y = r sin 0 sin <J> (38b) 

z = r cos 0 (38c) 


Since 


k = ka + k a + k a , 
— x— x y— y z— z 


then the "phase" term k * r in (19b) becomes 


k*r = kx + ky + kz 
— — x y J z 


= (k sin 0 cos <j> + k sin 0 sin <j> + k cos 0)r (39) 
x y z 


where k = J k 2 - k 2 - k 


Using the Euler identity. 


"l k * r 

we can also express the term e J — — in (19b) as 


e- 5 — — = cos (k • r) + j sin(k • r ) . 


The "phase" term k • r is in general composed of real and 
imaginary parts, since k , in (39), can be real or imaginary 
as seen from (20) . Thus letting k • r be described by 


jk * r = -y'r = -(a' + jg')r, 

\ 

where a' and B 1 are real numbers (attenuation and phase 
constants, respectively), we then get 

jk • r -a 'r -jg 'r 
e— — - e e . 


= e a r (cosB'r + j sing'r). 


(40) 
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When r is very large, we see from (40) that the real and 


imaginary components of e 


jk • r 


oscillate very rapidly 


between equal positive and negative values except for 
certain values of k and k for which k • r doesn't 

x y - 

change much for changes in k x and k^,. For those values 
of k^ and k the first order change in k • r is zero 

/N 

for small changes in k and k . Conversely A(k ,k ) 

^ x y — x y 

in (19b) is a slowly varying function in k and k (as 

x y 

compared with k • r being a rapidly varying function in k 

" X 

and ky) . Thus the contribution to the value of the integrals 

A ”1 )c * X* 

over all k x and k^ of A multiplied by e J — — comes 

A 

essentially only from the values of A determined by the k x 

and ky which make the ''phase”, k * r, "stationary" 

(constant) . The "phase" is stationary when 


3k • r 


3k 


3k 


= 0 


(41) 


Using the value for k * r from (39) , 


3k • r 
~Tk. 


■rz — (k sin 0 cos 4 + k sin 0 sin <f> 
3k x y 


+ / k 2 ._ , 2,2 


k - k cos 0 ) r = 0 
x y 


(sin 0 cos <j> — 


k„cos 0 


-) = 0 . 


- 23 - 


Therefore 


k sin 0 cos <p 

V = -1 

X cos 0 


(42) 


Likewise, 3k * r/3k yields 


k sin 0 sin d> 

k = _z 

y cos 0 


(43) 


But 


, 2 , 2,2 ,2 
k = k . - k - k , 
z x y 


which upon substitution of the values of k and k into 
^ x y 


(42) and (43) gives 


k 2 = k 2 - k 2 sin2 9 
Z z cos 2 0 


which yields 


k = k cos 0 . 
z 


(44) 


Combining the results of (44) with (42) and (43) yields 

k , k , k , which we will define as those value of 
x y z 
o 1 o o 

k x » ky, k z for which k • r is stationary, namely: 


k 

z 

o 

= k 

cos 

9 



(45a) 

k 

X 

o 

= k 

sin 

0 

cos 

0 

(45b) 

k 

V 

= k 

sin 

0 

sin 

<l> 

(45c) 
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(Note from (45a) that k = k is real and therefore 

z z 

o 

a = 0 in (21). Therefore we have non-evanescence.) The 
phase k * r can be expanded in a Taylor series about k 

x o 

and k : 


3k • r 


3k ‘r 


r = k * r 


k ,k 
x y„ 
o J o 


3k 


Ak + — — 
x 3k 

k , k 
x y 
o J o 


Ak 


k ,k 
x y 
o 1 o 


3 k • r 


+ 35 


3k 


2 1 3 - ‘ - 

Ak + ~ 

x 2.2 

i 1 k 

k ,k y 

x y 1 

o 1 o 


Ak 2 

y 

k ,k 

^ y. 


3 k 


3k 3k 
x y 


Ak k + 
x y 

k ,k 

x y 
- 1 


(46) 


where 


Ak = k - k 

XXX 


and 


Ak = k - k 

y y y 0 


But (41) tells us that 


3k • r 


3k. 


k ,k 


■ y 
O 1 o 


3k • r 
3k 


k ,k 

y 

o 1 o 


= .0, 


and (39) and (45) give us 
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k • r 


(k sin 0 cos <{> + k sin 6 sin <t> 


k ,k 
x y 
o 1 o 


+ /k 2 - k ^ - k ^ CQS 


x y 

O 1 o 


2 2 2 2 2 
= (k sin 0 cos cj> + k sin 0 sin 4> + k cos 0)r 


= kr 


Therefore, (46) becomes: 


k * r = kr + (aAk 2 + bAk Ak + cAk 2 ) (47) 

r~ x x y z 


where 


9 2 k * r 


a = h 


9k 


b = 


3 2 k * r 


k • r 

r C = h j , 

9k 


all evaluated at k = k and k = k 

x x o y Y o 

Therefore we have in (19b) that 


E(r) = E (x,y , z) — =■ 

4tt 


Atk ,k )e 

— 00 ^ —00 -* 


jk • r 


dk dk 
x y 


4 TT'' 


A (k ,k ) e 
-co- X o y o 


j (kr + aAk 2 + bAk Ak + cAk 2 ) 
J x x y y 


dk x dk y 


7T i (k x ' k y )ejkr 
4tt o-^o 


-00 J 


2 2 

j (au + buv + cu ), , ,. os 
e J dudv (48) 


where u = Ak and v = Ak . It can be shown that the double 
x y 

integral of the exponential term in (48) evaluates to 
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2 TT 


( 49 ) 




When r is very large, we are in the far field (subject 
to (1) ). In the far field then, we have (from (48) and 
(49) ) : 

~ jkr 

E(£>= j 2i — k 2 A< k x < k > < 5( » 

o o J o 

with k ,k ,k given by- (45) . Using (23) we can express 
o J o o 

A 

A in (50) in an explicit manner to give the far field 
in terms of the near field: 

. , jkr jkz cos0 

E(r) = j e o cos0 

,00 ,oo ^ 

x -jk sin 0(x cos cj> + y sin d) ) , , 

• E(x,y,z)e J v 1 T dxdy (51) 

J _oo — oo ® 

where, from (37) and (38) , 

r = a r sin 0 cos <b + a sin 0 sin d> + a r cos 0. 

— —x — y — z 


2.1 Probe Compensation 

Although (51) describes the far field precisely in terms 
of the near field, in a practical sense it does not 
describe the far field in terms of near field measurements . 
The reason of course is that any measurement of the near 
field, plane wave or otherwise, is altered by the measure- 
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ment probe characteristics. So the near field in a plane, 

A 

represented by E(x,y,z o ) in (51), would have to be 
obtained by deconvolution from actual near field measure- 
ments. Alternately, the far field can be obtained directly 
in a form that uses the near field measurements, but in 
conjunction with a description of the probe antenna pattern 
such that compensation is made for the probe characteristics. 
Such is the form given in References 3 and 12 and derived 
below. (However the derivation here is different than 
that in Reference 3 where the Lorentz reciprocity theorem 
is used instead.) 

Consider the geometry of the antenna and measurement system 
shown in Figure 4. For the moment, let the measurement 



x,y plane = antenna aperture plane x',y' plane = measurement plane 

= probe aperture plane 


FIGURE 4 

GEOMETRICAL RELATIONSHIP OF ANTENNA AND MEASUREMENT PROBE 
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system be an "infinitesimal" probe located at the coor- 
dinates x q , y q z ^ in the x,y,z coordinate system (0,0,0 
in the x' , y', z' system.) Let the probe's measurement 
characteristic be described by the vector/phasor term 

A A 

C = cdxdy, where 

A A A /S -A 

c = c(x,y,z) = c. (x , y , z) a + c 9 (x,y,z) a + c-, (x ,y , z) a . 

x x z y j z 

(52) 

When the probe is oriented so that it "points" in a 
direction perpendicular to the antenna aperture plane, 

A A 

we have C' = c'dx'dy', where 


c ' (x ' ,y ' , z ' ) = c x (x' ,y ' ,z ' ) a_ x , + c 2 (x',y',z') a y , 

+c 3 ( x ' , y ' ,z')a z , (53) 

But it can be shown from Figure 4 that the unit direction 
vectors and the differentials dx' , dy' in the x' , y', z' 
coordinate system are related to those in the x,y,z system 
as follows: 



dx' = -dx 
dy' = dy 
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Therefore a voltage received by the "inf itesimal " probe 
can be represented as 


dv( V y o'V = [-E x (x o , y o ,z o )c 1 (0,0,0) 

+ E y( X O 'y o ' Z o )C 2 (0,0,0) 

A A 

- E 2 (x o/ y o ,z o )c 3 {0,0,0)]dxdy 

A A 

= (VWV ^x + WW % 

+ WV z o>S Z ] • -[-Ci (0,0,0) a x 

A A 

+ c 0 (0,0,0) a - c., (0,0,0) a ] dxdy 
4 y o . z 

/V A 

= t E X (X 0’ y 0, Z 0>- + VWV 2y 

+ E z (x o' y o' z o , izJ ■ [^ 1 (0,0,0) a x . 

A A 

+ c 2 (0,0,0 ) a y ,+ c 3 ( 0 , 0 , 0 ) a z i] (-dx’)dy' 
= E(x Q/ y o ,z o ) • [-£’(0,0,0)] dx'dy'. (54) 
Equation (54) becomes 


dv(x o ,y o ,z Q ) = E(x,y,z) • [-c* (x 1 ,y ' ,z' )J dx'dy' 


where, from Figure 5, 


x'=y’ = z'=0 (55) 


x = x q - x' , y = y Q + y ' , z = z q - z’ . (56) 


Thus from (55) and (56) we get 


dv< VW 


t30- 


X 


E (x - 
— o 


A 

Y o + y '' z o _z ' 5 * [“£' ( x * »y ' > z ' )] 


dx '• dy 


X 1 =y ' =Z 1 =0 


( 57 ) 


Let a probe with an aperture area, Q , be now placed about 
the point 0,0,0 in the x',y' plane with the probe aperture 
parallel to the x,y plane, as shown in Figure 5. Let the. 


0, probe aperture, parallel 
to x - y plane 


ANTENNA 

FIGURE 5 

RELATIONSHIP OF PROBE APERTURE TO ANTENNA 

probe aperture be filled with an infinite number of 
infinitesimal probe antennas which are described by (57) . 
Let us further assume that each of the infinitesimal probe 
antennas responds equally to the same field (ie. , 
uniform illumination taper for the probe aperture). Then 

A 

from (57) the probe's response to the field E is 



v,x o’ y o' z o > ‘ 


dv 


E (x q - x ' , y q + y 1 / z Q - z ' ) • [— c ' (x ' ,y 1 , z 1 ) jdx'dy 
CKO, 0,0) 


z' =0 


= [ [ E (x - x', y + y', z ) * [-c ' (x ' ,y ' , o) ] dx ' dy ' (58) 

jJ — O Q O u J 

a 

If the probe aperture is not uniformly illuminated, then 
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v(x o' y o' z o ) 


E(x - x' , 
— o 


+ y'» 


z o> 


a 


* [-c ' (x 1 ,y ' ,o) f (x ' , y ' ) ] dx 1 dy ' (59) 


where f(x , / y') is the illumination taper. In (59), let 

A 

the term £ be defined such that 


£(x',y') = -c(x’ ,y' ,o) f (x' ,y*) . (60) 

A 

Thus we see that £ is proportional to the aperture field 
distribution of the probe. Now the illumination taper, 
f, is zero outside Q . Therefore (59) can be expressed 


as 


v(x , y z ) = 
o a o, o 




I (x 0 - x '- y o + y ’’ V 


• g ( x ' , y ' ) dx ' dy ' . 


(61) 


Let the probe aperture not be restricted to the location 

about x ,y . Let it be free to move about in the measure- 
o J o 

raent plane at z = z q . Therefore in general, the measure- 
ment probe response in the measurement plane is 


v(x,y, z .) 


E(x - x’ , y + y' , 


z o } 


£ ( x ' , y ' ) dx ' dy ' . (62) 


Let us examine the nature of the integrations performed 
in (62). To that end, let a variable, w, be the result 


of the convolution of two other variables, u and v, such 


that 

w(t) = f u(t - x)v(x)dx. 

J — CO 

Therefore the Fourier transform of w,^{w}, is given by 
S?{w} = &{u}9{v} . 

On the other hand let w be the result of the correlation 
of u and v such that 




Therefore the two-dimensional Fourier transform of the 

A 

voltage, v(x,y,z o ), measured by the probe in the near 
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field, where v is expressed by (62) , can be expressed in 

/\ 

terms of the two-dimensional Fourier transforms of E and 

2 as 



r°° - - j (k x x + k y) 

v (x ,y , z ) e y 

— 00 


dxdy 


E (x,y , z q ) e 


- j (k x + 
J x 


k,y) 


dxdy 


- -j (k x x + k y) 

cj_(x, -y)e y dxdy 


— 00 J — CO 


[AOc^kyJe 3 2 °J- G(k x , -k y ) 


(63) 


where the bracketed term is from (22) . The term G is 
equivalent to the far field radiation pattern of the 
probe. This is because the far field radiation pattern 
of an antenna results from the Fourier transform of 
the aperture field distribution which here is propor- 

/v 

tional to g. 


Multiplying both sides of (63) by the term 


e jkr 
2irr k z 

and evaluating the resulting expression at the modes given 
by (45) : 
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k = k = k sin 0 cos <j> 
x x_ 


k = k = k sin 0 sin d> 

y y 

J J o 


k = k - k cos 0 
z 


yields 


[j ! 


jkr 


jk z 

^ z o ~ 

„ k A (k ,k ) e ° ] • G(k ,-k ) 

2nr z — x y J — x y 

o o J o o J o 


jkr r°° 

* ihtz k z v(x,y ,2 ) 

O J —CO J — 00 


-j (k x x + k y y) 
e ° ° dxdy . ( 64 ) 


Upon substituting the result of (50) into (64) we get 


E (r) • G (k , -k ) 

— — — x y 

o ^o 

e jkr 3k z o Z o f°° r™ 

= => k z ~ 2 ^ e 
o 


-j (k X + k y) 

/ n o y o , 

v ( x , y , z ) e dxdy . 

-<x> J -oo ° (65) 


It is helpful to express this result in terms of spherical 
coordinates : 


E (r) 


A 

= E r (r , 0 ,<p) a r + E (r,e.,<f>)a 0 

far field 


+ E (j) (r,0,4>)a (j) . 

But the radial component of an electric field in the 

far field of an antenna is, for practical purposes, zero. 

Therefore, 
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E r (r,e,cf>) = 0 . 


The shape of the other components of the electric field, 

✓V /N 

E 0 and E^ are independent of the distance, r, to the 
antenna (by definition of the far field) . Thus without 
loss in generality we can write 


E (r) 


far field 


= E 0 (;0,<t>)a ft + E A (0,<j») a A . (66) 


-0 “<j> 


-<t> 


Likewise for G in spherical coordinates 


G(k fk ) 
x o y o 


= Q(k sin 0 cos <J> , - k sin 0 sin <fi) 

= G(k sin 0 cos (-<j>) , k sin 0 sin (— <f> ) ) - 
Vie see from Figure 5 that 


<j>' = -<)), 0« = e. 


(67) 


Therefore 

A A 

G(k , -k ) = G(k sin 0' cos cf> ' , k sin 0' sin <j) ' ). (68) 
x o y o 

A 

The far field radiation pattern, G, can then be rewritten 
as G(d’,<p') which can be expanded in spherical coordinates 
as 


GO' ,<*>•') = Gjte',*')^' + 6^(0' ,4.')^' 


+ G r (0 , ,* , )a r '.. 


(69) 
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Since the radial component of the electric field trans- 
mitted by an antenna is effectively zero in the far 
field, by reciprocity the radial component of an antenna's 
far field receiving pattern is zero. Therefore in (69) , 

A 

G r ( e ’ ,4>') = 0. (70) 

Adding this to the fact that from Figure 4 

V = - a e' v = v • (71) 


Equation (69) , with the use of (67) , becomes 


G - -G Q (0 ,-<p) a Q + G^ (8 ,-<)>) a^ 


Therefore the expression in (65) can be evaluated in 
terms of the angles 0 and <{> to give us the relationship 
between the far field pattern of an antenna and the far 
field radiation pattern of a measurement probe. The 
relationship is given in terms of the spatial Fourier 
transform of the voltage measured by the probe in the 

A 

antenna's near field and the vector components of E and 

A A A A 

”E 0 ( 0 , <J> ) G 0 ( 0 , -<f> ) + E^ ( 0 , (j) ) G^ (0 , -<J) ) 


= 3 


, jkr -jkz cos 8 
.ke J „ o q 

e cos 0 


2ur 


— 00 J — . CO 


V 


(x, y , z )e' jk sin 0(x cos + + y sin 9, dxdy. {72) 
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< oi 


The probe measurement will in practice be a set of discrete 
samples. The double integration in (72) then becomes a 
double summation discrete approximation which can be 
efficiently evaluated using the Fast Fourier Transform. 


Equation (72) expresses the far field radiation pattern 
of the antenna's electric field in terms of the two 

A A -A 

orthogonal field components, Eg, and E^. The terms G 0 

A 

and G(j) are the two known (orthogonal) components of the 
far field radiation of the probe. They are expressed in 
the direction (0, -<!>) from the probe to the phase center 
of the antenna (which is opposite to the direction (0,c(>) 
from the antenna to the probe.) Thus we have one (linear) 
equation and two unknowns. Therefore a second, linearly 
independent equation is needed to solve for the two 
unknowns (Eg and E^) . A rotation of the probe, which 
thereby changes the polarization axis (axes), provides 
the needed second equation. The second equation has the 
form 

A A A A 

- Eg(e,<j>)G Re (0,-c}>) + E (0,4>)G Re (6,-4>) 


ke jkr - jkz Q cos 9 
^ 2 nr 6 


cos 0 


— OO J —00 


; R (x,y,z„)e-5 k sin 0 (x cos ♦ + V sin ♦) dxdy (73) 


where the subscript, R, refers to quantities resulting 
from the probe rotation. 
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Therefore for each 0, <J> combination we have two independent 

A /*v 

equations for the two unknowns, Eg and : 


A A A A 

~ E 0 G 0 + V<D = 


E 0 G R0 + E 0 G R4> 


result of measurement #1 
= result of measurement #2. 


These two equations immediately above are solvable if 
the determinant 


V 


/\ /N 



g r0 G r4> 


g 0 g R4) 


g * g r0 * 0 


(The determinant will be zero for the case where the probe 
is circularly polarized because such a probe will obviously 
give no new information upon rotation.) 


2.8 Other Modal Expansions 

Our focus in the report has been the plane wave modal 
expansion - antenna patterns from measurements made over 
a planar surface. Far field antenna patterns can also 

9 

be obtained from measurements made over cylindrical 
and spherical surfaces"'' (Figure 6) . These surfaces 
arise as a natural consequence of the type of probe 
scanning employed, as does the planar. The cylindrical 
results from vertical scanning of the probe together 
with rotation is azimuth of the antenna under test. 
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CYLINDRICAL-WAVE MODAL EXPANSION 



SPHERICAL-WAVE MODAL EXPANSION 


FIGURE 6 

OTHER MODAL EXPANSION MEASUREMENT SURFACES 
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Spherical is employed by fixing the probe in space and 
rotating the antenna in azimuth and elevation. This 
results in simpler probe calibration. These scanning 
surfaces naturally give rise respectively to cylindrical- 
wave and spherical-wave modal expansions. The Fast 
Fourier Transform can still be used for the cylindrical- 
wave modal expansion, however it cannot be used for the 
spherical-wave 1 . The spherical-wave modal expansion is 
computationally the most complicated. 


3.0 Acoustic Waves 


Let us now direct our attention to acoustic waves. 
Consider the acoustic wave equation for the instantaneous 
pressure increment, p: 




a 

2 

where V is again the Laplacian, and c 

a 

sound in air (for the sonic simulator) . 
chromatic" (i.e., single tone) field, 


is the velocity of 
For a "mono- 


p = Real {pe^ wt } 


and therefore. 


V^p = 




= -k 


2 ~ 
P 


(74a) 

(74b) 
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where 



and X is the wavelength in air for a given radian frequency, 

ci 

a) = 2irf. In rectangular coordinates. 


/\ A 


p = p(x,y , z) , 


and, of course, 


2 2 2 
2 = + 

v 2 2 2 
3x z dy 3 z z 


Taking the two-dimensional Fourier transform of the wave 
equation (74) with respect to x and y yields: 


9 - -j(k x + k y) 

V^p(x,y,z)e y dxdy 


— OO ' —00 


= -k B(k ,k ,z) 

x y 


(75) 


where 


,v r™ ~ 

B(k x ,k ,z) = p(x,y,z) 

* — 00 ' —00 


-j(k x + k y) 


y dxdy . (76) 


Note the similarity between (75), (76) and (7), (8). We 
see that the development here is identical to that of the 

/s /s. 

electromagnetic case except that vectors such as E, B, 

A 

and A are replaced respectively with scalars such as 
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p, B, and A. Thus we obviously obtain the results below 
for the acoustic case. In the far field, then. 


p(r) = j 


. e 


jkr 

- k A {k , k ) 
2nr z x y 
o o ■ r o 


(77) 


where 


A (k ,k ) = e 
x y 
o 1 o 


- jk z 

J 2 O'' 00 f0 ° 


— 00 * —CO 


-j(k x X + k y y) 

p (x,y , z q ) e ° ° dxdy (78) 


and where 



= k sin 0 cos <p 
= k sin 0 sin <J> 


k = / k 2 - k 2 - k 2 = k cos 0 . 
z o x o ^o 

Probe compensation then yields 

jk z 

~ z o ~ 

A (k ,k ) e ° G (k , -k ) 
x ' y x y ' 

o 1 o o*o 


(79a) 


(79b) 


— CO ' — oo 


v ( x , y , z Q ) e 


-j (k x 
J x 

o 


V y) 


dxdy . (80) 


Multiplying both sides of (80) by the term 
e jkr 

j | k 

J 2nr z 

o 

and then rearranging terms gives 
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/V 


jkr 

■ e 

- 1 2trr 


e 


p(r)G(k ,-k ) 

x o y o 



z 

o 


k 


z 

o 


v (x , y , 


z 

o 


-j (.k 

) e 



+ k y) 

y 

° dxdy ( 81) 


which then becomes 


/s , jkr -jkz cos 0 

p ( 0 , <j)) G ( 0 , -<J>) = j- |— e ° cos 0 


v(x,y , z ) e - -^ sin 6<x cos<t,+ y sin ^ dxdy . (82) 


— CO' — CO 


Note that (82) is identical to (72) with the exception of 

A A A A 

course that p and G are scalars instead of vectors E and G. 
If we have chosen our coordinate system such that the z 
axis was parallel to one of the axes of polarization of the 
antenna instead of perpendicular as in Figures 2 and 3, 
we would get an equation in E 0 G 0 and EqGq similar to (72) 
and expressed in the resultant newly defined 0 and (j)'*'. 
However now the electric field in any plane where <)> is 
constant would be due only to that one particular polari- 
zation. The electric field perpendicular to that plane 
would be due only to the other orthogonal polarization 
component. This means that for linear polarization in a 
direction parallel to the new z axis, E^ would be zero 
in our new equation that is equivalent to (72). Since 
any electromagnetic wave can be decomposed into two 
orthogonal linear polarizations, the polarization component 
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perpendicular to the first would result in an - 
only component (Eg = 0) in our new equation equivalent 

/s 

to (72) . Thus our newly defined Eg and E^ are 
separable by polarization. Therefore acoustic simulation 
of electromagnetic waves .can be performed wherein the 
new 0 represents one polarization and the new <j> the other. 

A 

The sonically obtained voltage, v, in (82) simulates the 

A 

electromagnetically obtained voltage, v, in our new 
equation equivalent to (72) . (The simulation would be 
actually performed as two separate simulations, one for 
each polarization. The results of the two simulations 
would then be orthogonally superimposed to obtain the 
final result. ) 

4.0 Conclusions 

(1) We see that near field to far field transformations 
are valid to obtain far field antenna information from 
near field measurements. 

(2) Acoustic simulation of electromagnetic near field 
to far field transformations is valid. 

* 

(3) Polarization is easily accommodated acoustically 
in the simulation of electromagnetic waves (including 
near field to far field transformations.) 


*except for certain types of reflections external to the antenna. 
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SECTION II 


MULTIPLE BEAM FORMING AND STEERING 


PAGE 




1.0 


Introduction 


1.1 Background 

Simulation software has been developed for the Retro- 

directive Phase Control Sonic Simulator to form and steer 

14 

up to five sonic beams simultaneously and independently 
Two methods of beamforming and steering were investigated. 
One was retrodirective phase control and the other was 
a mathematically synthesized open loop phase control. 

1.2 Beamforming with Sonic Simulator 

Beamforming requires that the proper array amplitude and 
phase tapers be generated to form the beam (s) at the 
desired focal point (s). These tapers dictate the amplitude 
and phase of each array transmitting element’s output. 
Beamforming can be though of as a two-step process: 

Step #1. For each beam, determine for each transmitting 
element or group of elements the difference 
in two propagation phase delays. One phase 
delay is that between the beam focal point 
and an array reference point (near the center 
of the array.) The other delay is that between 
the beam focal point and the transmitting element 
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or group of transmitting elements. This 
process step can he called "phase tuning". 

Retrodirectiye Phase Control - In this mode, 
the Sonic Simulator's transmitting elements 
are grouped together four-each in subarrays. 

The propagation phase delay difference in 
the Sonic Simulator system is measured directly 
as the phase difference between the pilot beam 
signal as received at the reference receiver 
on the array and as received at the center 
of each subarray. The focal point of the beam 
subsequently transmitted by the array is 
determined by the focal point of the pilot 
beam source. 

Mathematically Synthesized Open Loop Phase 
Control - The phase taper is mathematically 
synthesized in this mode. The propagation phase 
delay difference is computed from the geo- 
metric difference in the path lengths, in 
wavelengths, to the beam focal point from the 
center of the array and from each transmitting 
element. It is easy to show that if, relative 
to the center of the array, the coordinates 

v 

of the beam focal point are desired to be at 
(X,Y,Z), and if those of a transmitting element 
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are at (x,y) , then the propagation phase delay 
difference is 

/ X 2 + y 2 + Z 2 = /(x - x) 2 + (Y . - y) 2 + Z 2 t A • 360°, 
where A is the wavelength. 

Step #2. Determine the phase and amplitude (or power) 
of the output of each array transmitting 
element. The multiple beam forming method 
used with the Sonic Simulator enables all 
transmitting elements to contribute to all 
beams . The phase taper and amplitude taper 
required to form each beam separately are 
first determined. Then the resulting phase 
taper and amplitude taper required to trans- 
mit all beams simultaneously are obtained by 
superposition — the phasor (vector) sum at 
each transmitting element of the phase and 
amplitude of its output which would be required 
to form each beam singly. This procedure is 
mechanized in the Sonic Simulator system for 
beams of equal power and shape by calculating 
each transmitting element's output phase and 
power (which is proportional to the square of 
its amplitude) according to Equations (58) and 
(59) and Appendix D. in . Reference 14: 
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N 

l sinA<j> 

A± ^ -1 n=l an 

A^c = “ tan ~ N ‘ 

I cosA* 
n=l 

9 N 9 

a|F(r) i [( l cosAcf) } + 

n=l an 


N 

( l sinAc(> 
n=l 

where 

A4> c = output phase of signal at transmitting 
element , 

N = number of beams. 

A^an = the P hases measured in Step #1 above. 

P = output power of signal at transmitting 
•N 

element . 

2 

a|F(r) | = power taper as a function of radial 

distance, r, from center of array. 

(Although the above phase and power equations 
were developed for the method of closed loop 
phase control^ 2 * * 5 , they are also valid for 
retrodirective and mathematical open loop 
phase control.) 

The amplitude (power) taper used in the experiments run 
with the simulator is a IQ dB four-step approximation 
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to a gaussian curve. The step sizes are 2.5 dB. This 
is the amount by which the simulator's array is designed 
to vary the transmitting element's output power levels. 

The program listing for the software used in the Sonic 
Simulator, including that for multiple beam forming, is 
given in 3.0 of Section II. The operation of the soft- 
ware is given in Reference 16. 

2.0 Experimental Results 

In the final analysis, the demonstration and evaluation 
of retrodireetive and mathematically synthesized open 
loop phase control techniques for multiple simultaneous 
beam generation with a phased array antenna is proven 
by the actual generation of those beams. This is easily 
done using the Retrodireetive Sonic Simulator to 
generate the beams and then making experimental inves- 
tigations with the Beam/Pattern Monitor and Display 
System. The resulting multibeam data that is taken 
with the monitor can then be used to evaluate each 
technique against standard antenna performance criteria 
and to evaluate each technique against the other to 
determine their relative performance. 

Figures 7 and 8 show far field patterns of multiple 
beams transmitted by the simulator array. The beams 
were focused in a plane 2.5 meters from the array. 
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FIGURE 8 


TRIPLE-BEAM FAR FIELD PATTERNS — MATHEMATICAL OPEN LOOP PHASE CONTROL 





("Left" and "right" in the figures are as viewed from 
the array.) The patterns are horizontal cuts made 
through the center of the beams in the focal plane. 

Each multibeam consists of three equal-power 10 dB 
truncated-gaussian amplitude taper beams. The top two 
beams were formed on centers about 1 m. apart and the 
bottom beam's center is about 0.8 m from the top two. 
The beams in Figure 7 were formed using retrodirective 
phase control and those in. Figure 8 using mathematical 
open loop. 

Although the sidelobes of the beams are higher than 
would be expected in free space, they are quite good 
and are to be expected for two reasons. For the sake 
of convenience, the measurements were made in a semi- 
anechoic environment. Secondly, the three beams in 
each multibeam were formed within several beamwidths 
of each other. The sidelobes of each individual beam 
are still of sufficient amplitude within that distance 
that at locations of phase reinforcement from two or 
more of the beams the resulting amplitude level can be 
much higher than that of a single sidelobe by itself. 

The mathematically formed beams are seen to be more 
sharply defined than the retrodirective ones with about 
3 dB higher amplitude and overall lower sidelobes. 

This is because for the sake of convenience the retro- 
directive beams were formed very close to the array 


- 56 - 


making it difficult to accurately illuminate the array 
with the retrodirective pilot signal. On the other 
hand, the mathematically formed beams were not as 
accurately pointed. A mathematically formed beam 
cannot track a receiver unless it has accurate infor- 
mation about the receiver location from another source. 

3.0 Source Code Listing - Simulation Software 

The following is the program listing, in BASIC, of 
the simulation software for the Retrodirective 
Sonic Simulator. 
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SECTION III 


RECOMMENDATIONS FOR ADDITIONAL AREAS OF WORK 


We recommend the following as additional areas of work 
to be performed by the contractor as part of future efforts 


1) Design and fabricate a scanner hardware and soft- 
ware control, data acquisition, and processing 
system to provide antenna patterns from measure- 
ment data output by the Beam/Pattern Monitor 

and Display System. A host computer, terminal, 
monitor, and printer would be supplied by the 
contractor as part of the system. 

2) Using the system developed in 1) above, develop 
a software system to generate far field antenna 
patterns from near fieid measurements made with 
the Beam/Pattern Monitojr and Display System. 

3) Provide maintenance and systems engineering 
support to NASA on the Sonic Simulator and Beam/ 
Pattern Monitor and Display Systems. 

4) Assist NASA in defining and investigating 
additonal Space Station antenna study areas. 

The effort could include but not be limited to 
multibeam phased array beamforming phase control 
using open loop, retrodirective , and/or inter- 
feromteric means. Such studies could also 
involve the use of the Sonic Simulator and 
Beam/Pattern Monitor Systems. 
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APPENDIX 


ENVIRONMENTAL EFFECTS ON SONIC SIMULATOR MEASUREMENTS 


A. 1 EFFECT OF TEMPERATURE ON SONIC SIMULATOR EXPERIMENTS 


The velocity of sound varies with temperature. Error 
can be introduced into near field measurements unless 
either the ambient temperature is held constant or 
else the frequency of operation is varied correspondingly 
with temperature to keep the wavelength constant. 

( The parameter of interest is phase - constant phase.) 
Therefore the sonic simulator is designed to vary its 
frequency output as a function of ambient temperature 
to maintain constant wavelength. 

The equation governing the relationship between the 
velocity of sound at any temperature is given by 

v = v / 1 + t/273 * (A-l) 

o 

where t is in °C, and v q is the velocity at 0°C. The 
wavelength. A, is related to the velocity, v, and 
frequency of sound, f, by the relationship 

A = -j- . (A-2) 


From the above two expressions we obtain the fact that 
at a fixed frequency, the wavelength changes about 0.2% 
for every 1°C change in temperature: 1 


Av/v ~ 1 . 1 

At 2 .273 


0.0018. 


♦Handbook of Chemistry and Physics, 23rd Edition, p. 1884. 
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For most far field measurements and for all of those 


beam characteristics measured as a part of previous 
contracts with the sonic simulator, that particular 
magnitude change in wavelength has negligible effects. 
For example, with an array of fixed dimensions, beam- 
width is proportional to wavelength. Thus, say, for 
a 10°C rise in temperature, there is only about a 2% 
increase in wavelength and hence the effect on measured 
beamwidth is negligible. Thus temperature effects on 
wavelength appear to have negligible consequences for 
far field measurements— those that normally are taken 
only of amplitude in the focal plane, a plane of 
constant phase. However, to determine far field pat- 
terns from near field measurements, phase must be 
preserved accurately. During the time it takes to 
make a complete scan of the near field, the phase error 
must be held to a fraction of a wavelength. This means 
that the change in wavelength due to temperature changes 
that occur within a complete scan of the sound field 
must be kept small. For example, consider measurements 
made in a plane parallel to the sonic array when the 
array is on the order of 30 to 60 cm. away from the 
plane. A measurement point which for instance is 60 cm. 
or so from transducers on the array is on the order of 

f\j 

7 wavelengths away from those transducers (1A = 8.6 cm. 
at 4kHz at room temperature.) For a 10°C change in 
temperature, we get a change of: 
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10 x 0.2% x 7A = 14% x A = 50° change in phase. 

Error can thus be introduced into near field measure- 
ments unless either the ambient temperature is held 
constant or else the frequency is varied correspondingly 
with temperature to keep the wavelength constant. 

Refer again to Equations (A-l) and (A-2) . If frequency, 
f, is varied according to the relationship 


f = f / 1 + t/273, (A- 3) 

then A will be constant and equal to v Q /f . Selecting 
f such that f = 4kHz at 25°C, and realizing that (A-3) 
is closely approximated by the linear expansion 


f 


V 1 + 



(A- 4) 


we see that if frequency is linearly compensated for 
temperature changes according to the expression in (A-4) , 
then wavelength is kept constant. The sonic simulator 
circuitry is designed to vary its output frequency in a 
linear fashion over a temperature range of approximately 
25°C (from about 10°C to 38°C) . 


A. 2 EFFECTS OF HUMIDITY AND BAROMETRIC PRESSURE 


The change in wavelength due to changes in humidity and 
barometric pressure is miniscule compared to changes due 
to temperature. This is because the velocity of sound 
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in air is such a weak function of those two variables, 
especially compared to temperature. 
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